arXiv:gr-qc/9808030vl 11 Aug 1998 


The Dirac particle on de Sitter background 


Ion I. Cotaescu 

The West University of Timi§oara, 

V. Parvan Ave. 4 , RO-1900 Timisoara 

February 4, 2008 


Abstract 

We show that the Dirac equation on de Sitter background can be 
analytically solved in a special static frame where the energy eigen- 
spinors can be expressed in terms of usual angular spinors known from 
special relativity, and a pair of radial wave functions. 


1 Introduction 


One of the most important but difficult problems of the quantum field 
theory in curved spacetime is the problem of finding analytic solutions of 
Dirac equation on given backgrounds. Since there are no general methods 
to obtain such solutions, these must be derived in each particular case sepa¬ 
rately, starting with a suitable chart and tetrad gauge and by using appro¬ 
priate calculation procedures. From the results reported up to now |1|. Q 
we understand that the form of the solutions obtained through separation of 
variables as well as their physical significance are strongly dependent on the 
choice of all these ingredients. In these conditions it is helpful to exploit, in 
addition, the effects of the global symmetries of the background. 

An important case is that of the Dirac equation in spherically symmet¬ 
ric (central) static charts which have the global symmetry of the group 
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T(l) ® SO (3), of time translations and rotations of the Cartesian space co¬ 
ordinates. Recently, we have proposed a Cartesian gauge 0 in which the 
whole theory remains covariant under this group when the Cartesian holo- 
nomic coordinates are used. Then the energy and angular momentum are 
conserved like in special relativity from which we can take over the method 
of separation of variables in spherical coordinates. In fact, our gauge defines 
Cartesian local (unholonomic) frames which play here the same role as the 
Cartesian natural frame of the Minkowski spacetime, since their third axes 
are just those of projections of the whole angular momentum. This allowed 
us to separate the spherical variables in terms of usual angular spinors such 
that all the constants involved in separation of variables get physical mean¬ 
ing. Therefore we are sure that the local properties of the Dirac held can be 
correctly interpreted. 

On this way we have found a complete formulation of the radial problem 
of the Dirac equation on central static backgrounds, obtaining the radial 
equations and the form of the radial scalar product in the most general case 
[[!]. Based on these results, we have solved the Dirac equation on anti-de 
Sitter static background, giving the formula of the discrete energy levels and 
the form of the energy eigenspinors in the spherical coordinates of a special 
natural frame. Moreover, we have observed therein that the solutions of 
the Dirac equation on de Sitter static backgrounds could be derived in the 
same manner, following step by step the technique of the anti-de Sitter case. 
This is just what we would like to present here. Our aim is to study the 
radial problem and its supersymmetry which helps us to find the radial wave 
functions we need in order to write down the energy eigenspinors in spherical 
coordinates (and natural units, h — c—1). 

We start in the second section with a short review of our main results con¬ 
cerning the separation of spherical variables in the Dirac equation on central 
static charts. The third section is devoted to the hidden supersymmetry of 
the radial problem which can be pointed out by using an appropriate trans¬ 
formation. This simplifies the form of the radial Hamiltonian allowing us to 
completely solve the radial problem, as shown in Sec.4. 

2 Preliminaries 
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The main point of our approach || is the choice of tetrads in Cartesian 
gauge and Cartesian coordinates such that Dirac equation takes a simple 
form which transforms manifestly covariant under the rotations of the space 
Cartesian coordinates. Consequently, its theory in the spherical natural co¬ 
ordinates, (t,r,9,4>), can be done like in special relativity. Our tetrad fields 
depend on three arbitrary functions of r, u, v and w, which give the line 
element 



( 1 ) 


and determine the form of the Dirac equation of the free field, ■?/, of mass M. 


We have shown that this equation has particular positive frequency solutions 
of energy E, 



( 2 ) 


Here UE,j,K :i ,m :i are the particle-like energy eigenspinors which depend on the 
radial wave functions and on the four-component angular spinors, , 

known from special relativity These are orthogonal to each other being 
completely determined by the quantum number j, of the angular momentum, 
the quantum number rrij, of its projection along the third axis of the local 
Cartesian frame, and the value of Kj = ±(j + 1/2). 

The radial wave functions are solutions of a pair of radial equations which 
can be written in compact form as the eigenvalues problem 


HE = EE 


(3) 


of the radial Hamiltonian 



H 


(4) 



in the space of the two-component vectors E = \f^ + \ /*■ / 7 where the radial 
scalar product is 
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This selects the ’’good” radial wave functions (i.e. square integrable functions 
or tempered distributions) which enter in the structure of the particle-like 
energy eigenspinors of (j2j). Furthermore, the negative frequency (antiparticle¬ 
like) eigenspinors can be obtained directly by using the charge conjugation 

i- 

Here we use three arbitrary functions (u, v, and w) which define the 
metric and implicitly the form of the radial Hamiltonian. However, it is 
known that, in general, for any central static metric two such functions are 
enough. This means that we have a supplementary degree of freedom allowing 
us to choose suitable radial coordinate. The most convenient is to work in 
the special natural frame £|] where, by definition, u — 1. According to ([[]) its 
radial coordinate is 

r s (r) — / + const. (6) 

J u(r) 

Here the constant assures the condition r s (0) = 0. 

3 The supersymmetry of the radial problem 


Let us consider now the problem of the massive Dirac particle on de Sitter 
background. There exists a static chart with usual spherical coordinates, 
(i, r, 0, 0), where the line element is 

ds 2 = (1 — a ) 2 f 2 )dt 2 -- — f 2 (d9 2 + sin 2 9 d0 2 ). (7) 

1 — to 2 r 2 


In the corresponding special frame, denoted fom now by (t, r, 6, 0), the radial 
coordinate is 


r = — arctanluuf, 

tO 

as it is obtained from (|) and ©• Here, the radial domain is D r 
and the line element has the form 

ds 2 =- ~ 2 — dt 2 — dr 2 -- sinh 2 tor ( d6 2 + sin 2 9 dtp 2 ) 

cosh' tor L to 2 


( 8 ) 
[0, oo) 

(9) 
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from which we can identify the functions v and w we need to write the radial 
Hamiltonian 

'■' h _ d i_ : 

Hr sinh< 


H = 


cosh 


our 


i our 


d 

Hr 


UJKj 

^ sinh our 


ouk 


( 10 ) 


cosh our 

with the notation k = M/ou (i.e. Me 2 /Tiuu in usual units). 

Our radial problem has hidden supersymmetry like in the anti-de Sitter 
case [fij . This can be easily pointed out with the help of the transformation 
T —> T = U{r)T where 


U(r) = 


cosh 
i sinh 


-i sinh ^rf- 
cosh 


( 11 ) 


A little calculation shows us that the transformed Hamiltonian, 


.ou 


H = U(r)HU~\r)-t--l 2x2 , 


which gives the new eigenvalue problem 


HH = E 


.ou 




has supersymmetry since it has the requested specihc form, 


H = 


d 

Hr 


_ d 
Hr 


+ W 


+ W 


-v 


Here v = ou[k — inj) is a constant and 


W(r) = ou(iktaahour + Kj cothcnr) 


( 12 ) 


(13) 


(14) 


(15) 


is the superpotential of the radial problem |4j. The transformed radial wave 
functions (which are the components of T) satisfy the second order 
equations resulted from the square of (|I3|). These are 


/ d 2 2 ik(ik ± 1) 2 K j( K j =t 1)' 
y dr 2 cosh 2 o;r sinh 2 our 

where we have denoted e = E/ou — i/2. 


/(±)( r ) = £n 2 e 2 / {± )(r), (16) 
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4 Solutions 


The solutions of Eqs.flTEl) are well-known to be expressed in terms of 
Gauss hypergeometric functions 0, F±(y) = F(a±, /3±,'f±,y), depending on 
the new variable y — — sinh 2 ujr , as 

f {± \y) = JV±(i - y) p± y s± F±(y) (17) 

where 

Z6 Z6 1 

T P± T T pz t ■^'7 7± 2s± T —, (18) 

N± are normalization factors while the parameters p± and s± are related 
with k and Kj through 

2s±(2s±-l) = K j (K j ± 1), (19) 

2p±(2p± — 1) = ik{ik ± 1). (20) 


Furthermore, we have to find the suitable values of these parameters such 
that the functions (|17D should be solutions of the transformed radial problem 
(p~3|) . If we replace ([L7|) in (|T3|), after a few manipulation, we obtain 


y (1 - y) ~ v (p± ± yj F ±(y) + (1 - y) (»± ± y) F±(.v) 

= ~ ± 1 ± ^- ±+1/2 (i - </) w -«+ i/ 2 K f (< / ), (21) 

where 77 = ±1. These equations are nothing else than the usual identities 
of hypergeometric functions if the values of s±, p± and IV + /IV_ are correctly 
matched. First we observe that the differences s+ — s_ and p + — p- must be 
half-integer since we work with analytic functions of y. This means that the 
allowed groups of solutions of ( |T9|) are 


2 S+ = ~Kj 
2 s__ — —Kj T 


1 , 


while Eq. ([201) gives us 
2 p\ 

2 p l _ 


—ik , 
—ik + 1 , 


2 Sy — Kj + 1, 

2s 2 = Kj, 


( 22 ) 
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2 p\ — ik + 1 , 
2p 2 _ = ik. 


(23) 





On the other hand, it is known that the hypergeometric functions of (|T7|) 
are analytical on the domain D y = (—oo,0], corresponding to D r , only if 
9 R( 7 ±) > 9fi(/3±) > 0 [| 6 ;]. Moreover, their factors must be regular on this 
domain including y = 0. Obviously, both these conditions are accomplished 
if we take 

s± > 0. (24) 

We specify that there are no restrictions upon the values of the parameters 
P±- 

We have hence all the possible combinations of parameter values giving 
the solutions of the second order equations (|T 6 | ) which satisfy the transformed 
radial problem. These solutions will be denoted by (a, b), a,b — 1,2, under¬ 
standing that the corresponding parameters are s±, p±, as given by (| 22 |) and 
), and 

^ n M „ h ( n ^ ~ 1 


(a,b) CL I \j i 

a± = s± + p ± + 


= s a ± +p b ± - 


7± } = 


+ 


2' -rt ■ jrzc 2 1 ,=c ’ 2 (25) 

The condition (|4[) requires to chose a = 1 when Kj = — j — 1/2, and a = 2 
if k — j + 1/2. Thus, for each given set ( E,j , Kj) we have a pair of different 
radial solutions, with b = 1,2. Therefore, it is convenient to denote the 
transformed radial wave functions ( 0 ) by bearing in mind that the 

value of a determines that of Kj. The last step is to calculate the values of 
N + /N_. From fl2lD it results 


N. 


(17) 


Kj — 


’ 3 2 


a 


(hi) 


IV. 


(1,2) 


Kj = j + 


IV. 

N. 


(1.1) 

f 

( 2 . 1 ) 


7?> 


V- 


p 


( 1 , 2 ) 


N. 


( 2 , 1 ) 


= 1 - 


0 1} 

( 2 ) 

7- 


IV. 

IV. 


( 1 , 2 ) 

4 - 

( 2 , 2 ) 


7+> 


a 


( 2 , 2 ) 


IV. 


( 2 , 2 ) 


( 2 ) 

7- 


1 , (26) 
(27) 


Notice that here 7 ^ = yi 2 " 1 — j + 1. 

Now we can restore the form of the original radial wave functions of (|2 
by using the inverse of ([11]). In our new notation these wave functions are 


SEjdS ) = COSll y/IjAiW ± ! 


(28) 


For very large r (when y —> — 00 ) the hypergeo metric functions behave as 
rs./ (— y) a± Thereby we deduce that f*E^ a ,b ~ exp(— ieur) and 

iEr 


f. 


(±) 

E,j,a,b 


(29) 
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This means that the functions (|28|) represent tempered distributions corre¬ 
sponding to a continuous energy spectrum. On the other hand, (^) indicates 
that this energy spectrum covers the whole real axis, as it seems to be natu¬ 
ral since the metric is not asymptotically flat. However, in our opinion, it is 
premature to draw definitive conclusions before to carefully study the prop¬ 
erties of the radial state space. Anyway, it is clear that the energy spectrum 
is continuous, without discrete part, while the energy levels are infinitely de¬ 
generated since there are no restrictions upon the values of j, which can be 
any positive half-integer. 

5 Comments 


Here we have derived the solutions of the Dirac equation on de Sitter 
background. This was possible grace of our general method based on the 
Cartesian tetrad gauge which preserves the maximal global symmetry of the 
central static charts. We must say that these solutions are different from the 
other ones obtained by using also Cartesian gauges but in the spherical co¬ 
ordinates of moving charts The argument is that our solutions rewritten 
in the comoving frame have no more separated variables. 

^From the technical point of view, we have used here the same procedure 
as in Ref. 0 . Thus it is clear that the Dirac equation on de-Sitter or anti-de 
Sitter backgrounds can be solved in the same manner. However, despite this 
fact, the results are different since the corresponding energy spectra are of 
different kinds. We believe that these are good examples of free Dirac fields 
with continuous and respectively discrete energy spectra which could help 
us to understand some sensitive aspects of the quantum theory in curved 
spacetime. We refer especially to the structure of the operator algebra and 
its dependence on the choice of the tetrad gauge in a given chart. A guide in 
this direction could be our recent study of the algebras of the one-dimensional 
relativistic oscillators []Tj] which have similar superpotentials as those of the 
transformed Hamiltonians of the de Sitter or anti-de Sitter radial problems. 
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